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In the experiment performed at PSI [1,2], isospin-symmetric S-wave vrA^ scat- 
tering lengths Oqj^ and a^^ are obtained^ from the measurements of the 
Ti~p atom strong energy-level shift eis{'K~p)., its decay width into n^n state 
ris(vr~p), and the strong energy- level shift of the ii'd atom eis{'K~d). These 
three measurements are claimed to be fairly consistent within the error bars, 
and - by using the technique described in Ref. [1] - yield the strong -kN scat- 
tering lengths in an accuracy which is unique for hadron physics: Oq"^ = 
(1.6 ± 1.3) X IQ-^M'l and a^^ = (86.8 ± 1.4) x 10~^M;^. The scattering 
length ao+ may be used as an input in the Goldberger-Miyazawa-Oehme [4] 
sum rule to determine the irNN coupling constant. Recently, a new experiment 
has been approved on the pionic hydrogen [5] that will allow one to measure 
^isiT^p) with a much higher accuracy and thus, in principle, to determine 



^ In Refs. [1,2], the quantities 69 = '^0+ ^^^ ^1 ~ ~'^o+ ^^^ used. Below, we follow 
the conventions of [3] 
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the TiN scattering lengths from the data on pionic hydrogen alone. This, in 
its turn, will enable one to vastly reduce the model-dependent uncertainties 
that come from the analysis of three-body problem [6]. Further, at CERN, the 
DIRAC collaboration conducts the experiment aimed at the measurement of 
the TT^TT^ atom lifetime within the 10 % accuracy that will allow one to deter- 
mine the difference of tttt scattering lengths a^ — 02 with a 5 % precision [7]. 
The DIRAC measurement provides a critical test for large/small condensate 
scenario in QCD: should it turn out that the quantity Oq — 02 is different from 
the value predicted in standard Chiral Perturbation Theory (ChPT) [8], one 
has to conclude that the chiral symmetry breaking in QCD proceeds differ- 
ently from the standard picture [9]. Finally, the DEAR collaboration [10] at 
the DA<I>NE facility (Frascati) plans to measure the energy level shift and 
lifetime of the Is state in K~p and K~d atoms - with considerably higher 
precision than in the previous experiment carried out at KEK [11] for K~p 
atoms. It is expected [10] that this will result in a precise determination of 
the / = 0, 1 S'-wave scattering lengths - although, of course, one will again be 
faced with the three-body problem already mentioned. It will be a challenge 
for theorists to extract from this new information on the ^A^ amplitude at 
threshold a more precise value of e.g. the kaon-nucleon sigma term and of the 
strangeness content of the nucleon [12]. 

The relations between the characteristics of hadronic atoms and hadronic 
scattering lengths in the leading order in isospin breaking are given by the 
well-known formulae by Deser et al. [13] - consequently, the measurement 
of these characteristics provides one with the unique possibility to directly 
extract these scattering lengths from the data. In order to perform the above- 
mentioned elaborate tests of QCD at low energy, the precision of the leading- 
order formula is, however, not sufficient. The theoretical calculations of the 
hadronic atom observables should be necessarily carried out at the accuracy 
that matches the accuracy of the experimental measurements. To this end, one 
may invoke the potential model, assuming, e.g. that the strong interactions 
which are described by energy-independent local potentials, do not violate 
the isospin symmetry. This framework, applied to the case of the tt+tt^ atom 
decay [14], failed, however, to agree with the magnitude and even the sign of 
the correction term to the leading-order formula, evaluated independently in 
ChPT [15-18]. The reason for this discrepancy is now well understood: the 
simple potential model described above, does not include the full content of 
the isospin-breaking effects in QCD. For this reason, in our opinion, it is legit- 
imate to question the applicability of the potential model to the calculation 
of the energy spectrum and decays of other hadronic atoms (e.g., the Ti~p 
atom), as well. These calculations should be also reexamined on the basis of 
low-energy effective theories where the isospin-breaking effects are completely 
taken into account. 

The present paper is aimed at establishing the precise relation between the 



strong energy-level shift of the TT~p atom in the Is state, and strong vrA^ scat- 
tering lengths. This relation, valid at next-to-leading order in isospin breaking, 
can be used to extract vrA^ scattering lengths from the measured energy spec- 
trum of the 7r~p atom. So far, only potential model calculations have been 
utilized for this purpose [1]. In the present paper, however, we use the non- 
relativistic effective Lagrangian approach and incorporate all isospin-breaking 
effects. The decay of 7r~p atom will be considered in a separate publication. 
Note that the same approach has been already applied to study the tt+tt" 
atom decay in ChPT [17,18]. 

Our strategy is the follows. First, we write down the general non-relativistic 
Lagrangian that describes the interactions between nucleons, pions and pho- 
tons at a very small momenta. This Lagrangian is then used - by means of 
Feshbach's formalism [19] - to calculate the ground-state energy of the tt~p 
atom. The energy is expressed in terms of the couplings of the non-relativistic 
Lagrangian which, in its turn, through the matching condition are related to 
the relativistic iT~p elastic scattering amplitude at threshold, calculated in 
ChPT. At the last step, we analyze the isospin-symmetry breaking contribu- 
tions to this amplitude at 0{p^) in ChPT that, using the already established 
relation between the amplitude and the bound-state energy, provides us with 
the isospin-breaking correction to the Is energy level of the 7[~p atom. 

The preliminary remarks are in order. The isospin breaking it the vrA^ sys- 
tem is due to two physically distinct sources: electromagnetic effects being 
proportional to the fine structure constant a, and up and down quark mass 
difference m^; — 171^. It is convenient to introduce the common counting for 
these two effects. In particular, one may define^ the small parameter 6 so 
that a ~ {rrid — rriu) ~ S. With this definition, both AM^ = 0{6) and 
nip — rrin = 0{S). Further, in the leading order in the fine-structure constant, 
the binding energy of the ir'p atom ground state is given by the well-known 
non-relativistic Coulomb-Schrodinger value Ei = — |yUca^, where fic denotes 
the reduced mass of n~p system. The leading-order strong energy-level shift is 
given by the strong amplitude times the square of the Coulomb wave function 
at the origin and, consequently, is of order 0{a^). We are interested in the 
leading order isospin-breaking corrections to the strong shift - thus we have 
to evaluate the bound-state energy at 0(6'^) that stands either for 0{a'^), or 
for 0{a^{md - m„))Q. 



^ Note that in Refs. [17,18], a slightly different convention a ~ (m^ — m„)^ ~ 5 
was adopted. The counting rule is purely a matter of convention: if we assume the 
same rule as in the present paper, the 0{{md — mu)'^) contributions which emerge in 
Refs. [17,18] , wiU be merely relegated to higher orders in the perturbation expansion. 
^ There is one exception from this counting rule. We include the effect due to the 
electron vacuum polarization as well. Formally this contribution is of order a^, but 
it is amplified by the large factor {M^+ /nief'- 



We emphasize that in this paper we deal with the shift of the ground-state level 
while the quantity that is experimentally measured is the transition energy 
between 3p — Is levels. The strong shift in the 3p state is, however, very small 
(see [1] and references therein), and for the time being is neglected in our 
approach. 

Now, we start with writing down the non-relativistic Lagrangian which will 
be used to calculate the energy of the n~p atom. This Lagrangian, in general, 
consists of an infinite tower of operators with an increasing mass dimension - 
it contains all possible terms that are allowed by discrete P, C, T symmetries, 
rotational invariance and gauge invariance. The convenient building blocks for 
this Lagrangian, therefore, are: the non-relativistic pion and nucleon fields, 
their covariant derivatives, and electric (E) and magnetic (B) fields. For the 
general procedure of the derivation, we refer to Ref. [20] which deals with the 
same problem in detail, in the context of QED, and merely display the result 
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-I- di ip^i/j ttItt- + d2{ip\ vtItto + h.c.) + d^x\ t^qT^o H • (1) 

Here F^^ stands for the electromagnetic field strength tensor, ip, x, -k^ and 
ttq denote the non-relativistic field operators for proton, neutron, charged and 
neutral pion fields, and 'Dtip = dfip — ieAQip, Vip = V ip + ieAip, Dt7i± = 
dtTi± =F ^e^o^i) D vT-i- = V TT-t ± ieATT-i- are the covariant derivatives acting 
on the proton and charged pion fields, respectively. The ellipses stand for 
the higher- dimensional operators and the counterterms that are needed to 
render the off-shell Green functions finite in the perturbation theory. The 



couplings Cp, Cp, Cp, c^, c^, cf , c^, di, d2, d^,- ■ ■ should be determined 
from matching of non-relativistic and relativistic amplitudes (see below). In 
general, the matching procedure determines the same combinations of the 
non-relativistic couplings that are needed in the bound-state calculations [21]. 

The Lagrangian displayed in Eq. (1), contains an infinite tower of operators 
with the increasing mass dimension. This Lagrangian is used in the perturba- 
tion theory to calculate the scattering amplitudes in the one-nucleon, one-pion 
and one-nucleon-one-pion sectors, with any number of photons, to a given ac- 
curacy in the coupling constant e and in the expansion in external 3-momenta. 
At any order in the perturbative expansion, these amplitudes can be renor- 
malized by adding the appropriate counterterms to the Lagrangian (1). More- 
over, these amplitudes coincide with the same amplitudes calculated in the 
relativistic theory, provided the couplings in the non-relativistic Lagrangian 
are matched in a sufficient accuracy in e. The perturbative calculations in the 
non-relativistic theory are carried out using the "threshold expansion" [22] 
in all Feynman integrals - this ensures that the counting rules are obeyed so 
that the higher-dimensional operators do not contribute to the momentum ex- 
pansion of the amplitudes at lower orders. In particular, this ensures that the 
mass parameters in the Lagrangian (1) coincide with the physical masses of the 
particles - there is no mass renormalization. Further, the same non-relativistic 
Lagrangian is used to calculate the bound-state characteristics. Here, the pow- 
ers of external momenta "translate" into the powers of e and the threshold 
expansion should be again used to ensure the validity of counting rules. For 
the details, we refer the reader to Ref. [21] and references therein. In actual 
calculations, one always deals with the truncated Lagrangian. In particular, 
in the present paper we aim at the evaluation of the ground-state energy of 
the 7r~p atom at the accuracy 0((5^). It can be demonstrated that then it is 
necessary to construct the non-relativistic Lagrangian that describes correctly 
all contributions to the ir^p elastic scattering amplitude at 0{6), which do not 
vanish at physical threshold. The Lagrangian displayed in Eq. (1), is sufficient 
for this purpose. It can be shown that all operators that can be constructed 
(including the counterterms), except the ones given in Eq. (1), will not con- 
tribute to the 7r~p scattering amphtude at 0{S) at threshold neither at tree 
level, nor through loops. 



Below, we consider briefly the non-relativistic couplings c^ , Cp, Cp, c^, c^, cf , 
c^. At the accuracy we are working, it is sufficient to match these constants 
at 0(e°). To this end, one may consider the nucleon and pion electromagnetic 
form-factors in the external field A^. The necessary details, including different 
normalization of states in the relativistic and non-relativistic theories, are 
provided in [20], where the same problem is treated in the context of QED. 
Below, we merely display the result of the matching 
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4 = 1 + 2/ip 

c^ = M^ < r^ > , (2) 

where /ip and /i„ denote the anomalous magnetic momenta of proton and 
neutron, respectively, and < r^ >, < r^ >, < r^ > stand for the charge 
radii of proton, neutron and charged pion. The remaining constants di, i = 
1 ■ ■ ■ 3 are determined from matching the ttN scattering amplitudes in different 
channels. 

The Lagrangian (1) still contains the terms that do not contribute to the 
ground-state energy level shift of the 7i~p atom in the accuracy we are working. 
In order to simplify this Lagrangian, in the Coulomb gauge we exclude the Aq 
field by using the equations of motion, and neglect high- dimensional operators 
that arise in a result of this operation. Further, we neglect the operators that 
contribute to the ttN scattering sectors with non-zero total charge, and to the 
spin-flip part of the 7i~p elastic scattering amplitude (since the ground state of 
the 7i~p atom is an S'-wave state, the spin-flip part does not contribute due to 
the rotational invariance). In a result, we arrive at the simplified Lagrangian 
which is better suited for the bound-state calculations 
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For the calculation of the energy-level shift of n'p atom ground state we again, 
as in Ref. [17] use the Feshbach's formalism [19]. To this end, we first obtain 
the total Hamiltonian of the system from the Lagrangian (3) by using the 
canonical formalism 



H = Ho + He + Hr + eH^ + Hs + e^AH^ = Hq + He + V , (5) 

where Hq is the free Hamiltonian describing photons and non-relativistic pions 
and nucleons. Further, Hr = / rf^xiJr, T = C, R,7, S, A, and 
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ifs = -c/i^t^7]-l7r_ -c?2(V^^X7rl7ro + h.c.) - d^X^xA'^o , 

Hx = i'H'nl'n-. (6) 

We treat the problem in the perturbation theory: Hq + He is the unperturbed 
Hamiltonian, whereas V is considered as a perturbation. The ground-state 
solution of the unperturbed Schrodinger equation in the CM frame 
(^1 - Ho - He)|^o(0; s)) = 0, with Ei = rup + M^+ + Ei, is given by 



l*o(0; ^)) = / (^ ^o(p) blip, s) al(-p) |0) , (7) 

where a_(p) and b\_{p,s) denote creation operators for non-relativistic n~ 
and proton acting on the Fock space vacuum, s is the projection of the proton 
(atom) spin, and ^E'o(p) stands for the non-relativistic Coulomb wave function 
in the momentum space. 

Now, we are going to evaluate the energy-level shift of the ground state due to 
the perturbation Hamiltonian V. To this end, we again, as in Ref. [17], look for 
the poles of the scattering operator T{z) = (He + V) + (He + \')Go{z)T{z) 
on the second Riemann sheet in the complex z-plane: the real and imagi- 
nary parts of the pole position give, by definition, the energy and the decay 
width of the metastable bound state. The free and Coulomb Green opera- 
tors are defined as Go(-2) = {z — Ho)^"*^ and G{z) = {z — Hq — He)~\ re- 
spectively. Further, we define the "Coulomb-pole removed" Green function as 
G{z) = G(z)(l — n), where H denotes the projector onto the Coulomb ground 



state \E'o (7). The ttN scattering states in the sector with the total charge are 
defined as |P, p; s)+ = fc^Pi, s) aL(p2) |0) and |P, p; s)o = 6o(pi, s) aj(p2) |0) 
(ao(p) and 6o(p, s) denote the creation operators for the tt" meson and neu- 
tron, respectively). The center-of mass and relative momenta are defined by 
P = pi + P2, P = ^2^^Pi - ^i^^P2, A = +,0, and r][^^ = nip/ {nip + M^+), 
77^+^ = M^+/{nip + M,+ ), r]f^ = nin/{nin + M^o), r^f = M,o/(m„ + M^o). We 
remove the CM momenta from the matrix elements of any operator R(-2) by 
introducing the notation 

a(P, q; s\R{z)\Q, p; s% = {2'af6\V) (q, s|r^B(z)|p; s') . (8) 

In general, the expression for the matrix elements of tab{z) contains the spin- 
nonfiip ad the spin-fiip parts that are defined in the following manner 

(q; s\yab{z)\^; s') = 5,,/(q|r^5(2;)|p) + i{cTss' ■ [p x q]) (q|r^^(2;)|p) . (9) 

Below, we define the "Coulomb-pole removed" transition operator that satis- 
fies the equation 

M{z)=Y + VG{z)M{z). (10) 

With the use of the Feshbach's formalism and the definitions given above, it 
is straightforward to demonstrate that the scattering operator T{z) develops 
the pole at z = z where z is the solution of the following equation 

z-E,-{^o\nil4z)\^o) = 0, (11) 

where (p|^E'o) = ^o(p) and m^^{z) is related to M(2;) through the defini- 
tions (8) and (9). 

In order to get the shift of the ground-state energy, the quantity m^_^_{z) 
is determined perturbatively from Eq. (10). The iteration series that emerge 
from Eq. (10), can be truncated since the higher-order terms do not contribute 
to the energy at 0(6'^). Using the explicit expression of V given by Eq. (5), 
replacing G{z) by Go(-2) whenever possible, and retaining only those terms 
that contribute at the accuracy we are working, the operator M(z) can be 
written in the form Nl{z) = \J{z) + W{z), where 

V{z) = Hr + e2H^Go(^)H^ + c^AHa , 
W{z) = Hs + HsG(z)Hs + HsGo(^)HsGo(^)Hs , (12) 



At the accuracy 0((5^), the energy of the bound state is equal to z = Ei + 
A^P + ei„ where 

AEr = Re (^oK+(^i)|^o) + E"^ , ei, = Re {^o\K+m\^o) ■ (13) 

Here u"^(2;), w"^_(2;) are related to U{z), W^z) through the definitions (8) 
and (9), and E^^'^ stands for the contribution due to the electron vacuum 
polarization (the corresponding term is added "by hand", see below). 

Using explicit expressions (6), we obtain 



Re(*„K,(^0l*„)^-f«VJ^^-;;5|;+4„V?A, (14) 

Re (vp„|w:;+(Bi)l*o) = — (-di + ^f ( « + — (l"" " 1)) + XT ''j* 

where /xq is the reduced mass of the ir^n system, and 



go = [2/io(mp + M^+ -m^- M^o)f^ (15) 



Here d and /x denote the dimension of space and the scale of the dimensional 
regularization, respectively (we use the dimensional regularization scheme 
throughout the paper). 

The vacuum polarization contribution, which is calculated separately, should 
be added to the electromagnetic energy shift. The modification of Coulomb 
potential due to electron vacuum polarization is given by the well-known ex- 
pression [16] 
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4»n| 

where irte denotes the electron mass. Applying the first-order perturbation 
theory, we find 



^vac ^ f dh^l{v)AV''^\r) = -^J^r]^<!>{r]), r] = ^^ , (17) 

J 671 rrie 

^ *(^) = ^(2vr-4r/ + -r/2 - -^'^j + ^,^^,_-^ Hv + ^V' " 1)- 
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Table 1. Contributions to the electromagnetic binding energy of the tt~p atom 
(eV). Higher-order corrections: vacuum polarization correction (order a^) and vertex 
correction have not been calculated. 

Corrections Notation Ref. [1] Our 

Point Coulomb, KG equation 

Finite size effect (proton, pion) 

Vacuum polarization, order a^ 

Relativistic recoil, proton spin and 

anomalous magnetic moment 

Vacuum polarization, order q^ 

Vertex correction 

The calculation of the electromagnetic energy-level shift is now complete. In 
order to compare our results with those given in Ref. [1], it is convenient to 
define 

El + A^r = E^^ + E^" + ^™' + ^'"^ , (18) 

where 



^KG 


-3235.156 


-3235.156 


_gfin 


0.102 


0.100 


^vac 


-3.246 


-3.241 


_g;rel 


0.047 

-0.018 

0.007 


0.047 



i^-^4.e«^(l + ^), 



i^'" = ^«V^((r^) + (r^)), (19) 

2^' U(mp + M^+) {mp + M^+y {nip + M^+] 

and E'™^ is given by Eq. (17). In Table 1 we collect various contributions to 
the electromagnetic energy-level shift calculated using the same input param- 
eters, as in Ref. [1]. As can be seen from Table 1, the calculated values of the 
electromagnetic shift almost coincide. 

In order to complete the calculation of the strong energy-level shift, one has to 
match at the accuracy 0{6) the particular combination of the non- relativistic 
couplings di, i = 1- ■ -3 that appears in Eq. (14). To this end, we consider the 
scattering operator 



Tr(z)=Vr + VrGr(z)Tr(z) 
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Vr = He + eH^ + Hs + c'AHa , Gr(^) = (^ - Ho - Ur)-' . (20) 

In the scattering operator Tr(z), all kinematical insertions contained in Hr 
are summed up in the external lines (see [21] for the details). We calculate the 
matrix element of the scattering operator Tr(z) between the 7r~p states at 
0{6). After removing the CM momentum, the spin-nonfiip part of this matrix 
element on energy shell is equal to 



fait" (Z)\U) - ^""" - ^^"^ (q + P)' , e2iaecm .„|Z„ . ^1 ^ 

(q|ti,,+4z)|p)- |^_p|^ 4m,M.+ |q-p|2+" (q|tR,++(^)|p) , 

(21) 

where the (divergent) Coulomb phase in the dimensional regularization scheme 
is given by 



^c(|p|) = f^/-^(j^-i(r(l) + ln47r)) + flln^, (22) 

IpI \d — 3 2 J \p\ n 

and 



^ / ,-r> / M N naiicdi anj.d'^ , Ipl o. 

Re(q|tJ^_++(^)|p) = ^ + _^^^ ini^ + e^A 

IPI ^ l^c 

-d, + ^dld, + dl^ + ---, (23) 

47r^ 

where ellipses stand for the terms that vanish at threshold. 

In order to carry out the matching, we consider the elastic 7r~p scattering 
amplitude T^^-p^n-p, calculated in the relativistic theory at 0(5)[^. This am- 
plitude can be uniquely decomposed into the piece containing all diagrams 
that are made disconnected by cutting one photon line, and the remainder: 

Tjr-p^TT-p = Tex + Ti, with 



Te. = J u{p^s) [l,F,{t) + ia,,{p, - q^Y ^ } u{q,s') {p2 + ga)^ /(t) , (24) 

where {qi,s), q2 and (pi,s'), p2 are the four- momenta and spin of outgo- 
ing/incoming particles, t = (gi —piY, and Fi(t) and f(t) denote, respectively, 
the nucleon and pion electromagnetic form factors. Then, from the matching 



'^ Note that, in order to have the same sign as in the non-relativistic definition, the 
sign of the relativistic amplitude is defined by 5 = 1 — iT. 
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condition T^-p_,^-p = (2Eg2Ep2wg2wpy^'^{q; s\tYi^^+{z)\p; s'), where Eg,Ep 
and Wq, Wp denote the relativistic energies of the outgoing/incoming nucleons 
and pions, respectively, one may conclude, that the spin-nonflip part T" of 
the remainder amplitude Ti should have the following threshold behavior 



Re I e-2-ec(|p|) rpA = ^ + ^2 In — - 87r(mp + M^+)A+--- , (25) 

L J I PI fie 

where p denotes the relative momentum of n^p pair in the CM frame, Bi, B2 
do not depend on p, and 



- ^^ = -rfi + ^ did, + dli (26) 

where A is the regular part of 7r~p scattering amplitude at threshold. 
Substituting Eq. (26) into Eq. (14), we finally arrive at 

ei, = -2a Vc ^ (1 - 2a/ic (In a - 1) ^) + ■ • • , (27) 

where the UV divergence contained in the quantity ^, has been cancelled. 

The equation (27) is the main result of the present paper - it gives the relation 
between the measured strong energy-level shift of the 7r~p atom and the 7r~p 
scattering amplitude at threshold, defined by Eqs. (24) and (25). This rela- 
tion is valid at 0(5), and in all orders in the chiral expansion. The threshold 
amplitude still contains the isospin-breaking effects. Consequently, in order to 
extract from the experimental data the vrA^ scattering lengths which are de- 
fined for the case with no isospin symmetry violation, one has to single out the 
isospin-breaking effects in the amplitude. To this end, one may invoke ChPT 
- below, we present the results of the calculations at the tree level. 

In order to evaluate the vrA^ scattering amplitude, we use the effective chi- 
ral pion-nucleon Lagrangian at O(p^) [23-25]. The amplitude is given by the 
following expression 



A = at^ + a^^ + e = ——^—-—-\mpM^ ^V^ 



+ mp(-8ciM2o + 4(c2 + C3)M2+ - Ae^ h - 6^/2) j , (2J 



where Cj, /, denote the low-energy constants in the 0{p'^) chiral pion-nucleon 
Lagrangian, qa is the nucleon axial- vector coupling constant [23-25], and in 
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our calculations we use the value F^r = 92.4 MeV. The scattering lengths 
^0+? '^0+ ^^6 calculated in the isospin symmetric theory with e = 0, m„ = ma 
and where, by convention, the masses of the pions and nucleons are taken 
equal to the physical masses of the charged pion and proton, respectively. With 
this convention, the isospin-breaking correction to the 7r~p elastic scattering 
amplitude at O(p^) is equal to 



^Sci(M^^-M^2o)-4eVi-eV2^ (29) 



7r{mp + M^+)F^ 



Below, we present the results of the numerical analysis for the isospin-breaking 
correction to the tt~p atom energy- level shift. In order to evaluate the second 
term in the brackets in Eq. (27), one may safely approximate A = aj^. + a^j^ 
and use the scattering lengths given in Ref. [2]: the total correction coming 
from this term amounts up to +0.66 x 10^^. Further, in order to evaluate the 
isospin-breaking correction to the n^p scattering amplitude at threshold given 
by Eq. (29), one has to specify the values of the O(p^) low-energy constants 
ci, /i and /2. We use the value of ci determined from the fit of the elastic nN 
scattering amplitude at threshold to KA86 data [26]: ci = —0.925 GeV~^ The 
value of the constant /2 can be extracted from the proton-neutron electromag- 
netic mass difference [27]: e^/2 = (— 0.76±0.3) MeV. The determination of the 
constant /i from data is however, problematic. For this reason, in our analysis 
we have used order-of-magnitude estimate for this constant: — I/2I < /i < |/2|- 

The prediction for the 7T~p ground-state energy-level shift based on 0{p'^) 
calculations in ChPT, strongly deviate from the predictions of the potential 
model. The prediction for the total isospin-breaking correction in the strong 
energy-level shift, defined according to [1] eis = eff^(l + 6^), where eff stands 
for the leading-order shift given by Deser formula [13], is: 6^ = (— 4.8±2.0)-10~^ 
at 0{p'^) in ChPT. This result differs more than twice from the prediction 
based on the potential model [1]: 5^ = (—2.1 ± 0.5) ■ 10~^, and the uncertainty 
caused mainly by a poor knowledge of the constant /i, is four times larger than 
the estimate of the systematic error in the potential model. In demonstration 
of the above discussion, in Fig. 1 we confront the results of the present analysis 
with those of the potential model, using the experimental data on hadronic 
atoms taken from Ref. [2]. 

We would like to emphasize that the physical effects encoded in the low- 
energy constants Ci, f\ and /2, are absent in the potential model. Namely, 
/i and /2 contain the effect of the direct quark-photon interaction, and the 
correction proportional to ci contains the effect coming from the dependence 
of the scattering amplitude on the quark mass - similar effect in the tt+tt" 
case constitutes the tree-level "mass-shift" correction (see, e.g. [15,16,18]). For 
this reason, the discrepancy between the potential model predictions and the 
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Fig. 1. Determination ofi^N scattering lengths from the pionic hydrogen and pionic 
deuterium measurements. Solid line corresponds to the energy-level shift calculations 
at 0{p'^) in ChPT, and the dashed lines - to the potential model resultsfl] 

results of the present analysis based on ChPT, does not come to our surprise. 
It remains to be seen, how the 0{p^) results in ChPT are altered by the loop 
corrections. A reliable estimate of the constant /i is also desirable. 

To conclude, we have derived the general formula that relates the energy-level 
shift of the vr^p atom to the vrA^ scattering amplitude at threshold. Numerical 
analysis carried out on the basis of this formula at O(p^) in ChPT, already 
indicates at the necessity to critically reaccess the values of the tiN scattering 
lengths, extracted from the energy-level shift measurement by means of the 
potential model-based theoretical studies. 
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